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Abstract 

We study minimum energy configurations of ./V particles in M 3 of charge 
— 1 ('electrons') in the potential of M particles of charges Z a > ('atomic 
nuclei'). In a suitable large-N limit, we determine the asymptotic electron 
distribution explicitly, showing in particular that the number of electrons 
surrounding each nucleus is asymptotic to the nuclear charge ("screening"). 
The proof proceeds by establishing, via Gamma-convergence, a coarse-grained 
variational principle for the limit distribution, which can be solved explicitly. 



1 Introduction 

The goal of this paper is to shed new light on basic screening effects in molecules. By 
screening one means the remarkable tendency of electrons to usually group them- 
selves around the atomic nuclei in such a way so as to cancel much of the long range 
~ 1/R Coulomb potential exerted by the nuclei and make the net potential exerted 
by the atoms short-range. 

Screening is usually tacitly assumed in molecular mechanics, molecular dynamics, 
statistical mechanics, and continuum mechanics. One starts from the outset from 
short-range atomistic forces respectively short-range continuum forces (i.e. stresses 
alias surface forces). 

Large-scale failure of screening (which is not observed in nature) would lead to 
spectacular breakdown of these models, e.g. bare Coulomb interactions violate the 
linear scaling of energy E with volume V for a quantum mechanical crystal. To see 
this, place bare atomic nuclei on the integer lattice points in a 3D cube of sidelength 
L, {R\, ...,Rm} = fl [0,L] 3 and evaluate their interaction energy asymptotically 
in the limit of large L: 

E ~ £ Ha,«^ \r^\^ !!\ xm < L \^[d{x-y)d{x+y) 

\Ra\,\R \<L 

~ L 3 -L 2 ~ M 5 / 3 ~ (1) 
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i.e. the energy per atom tends to infinity as the system gets large (in a finite system 
of M = 10 23 atoms it is already too large by a factor of about 10 17 ). 

The above example is unstable, but small-scale failure of screening is common in 
nature, and yields important 0(1) contributions to the energy per atom. Examples 
include ionic crystals like NaCl, molecules with low permanent multipole moment 
like H 2 0, intermediate states during chemical reactions, and core regions of atoms. 

We know of no mathematical results which directly explain and quantify screen- 
ing from full quantum mechanics. The perhaps furthest result in this direction 
concerns an indirect, coarse-scale manifestation of screening: the ground state en- 
ergy of a molecule with M atoms is known not to scale like the example (pQ), but is 
bounded above and below by a constant times M [DL67, LD68, LT75]. 

More insight has been obtained in asymptotic limits. For atoms in the limit 
of large atomic number Z, it is known [LS77] that the total electron density is 
asymptotically radial, the profile given by Thomas-Fermi theory, and falls off like 
r -6 , and so the asymptotic net potential exerted by the atom is short-range. For 
a closely related result see [ILS96]. Another interesting limit is the thermodynamic 
limit for crystalline solids, in which the nuclei are arranged on a regular subset 
of a periodic crystal lattice, say Br fl Z 3 , where Br denotes the ball of radius 
R around the origin, and R tends to infinity. In this case, for a slightly simplified 
version of quantum mechanics (absence of spin, a rigid wall assumption, and coupling 
to an electron reservoir), the ground state energy is known to be asymptotically 
proportional to the number of nuclei [Fe85]. Moreover for convex density functional 
models such as the Thomas-Fermi- Weizsacker model, the ground state density is 
known to become asymptotically periodic [CLL98]. The latter result, by Catto, Le 
Bris and Lions, can be viewed as a quantitative version of screening: the asymptotic 
amount of electron density in each unit cell exactly cancels the nuclear charge in 
that cell, making the net electrostatic potential excerted by the cell short-range. 

Here we introduce and analyze a model which allows some new insight into 
screening for general, non-periodic, arrangements of nuclei, at the expense of further 
simplification of the treatment of electrons. The model maintains the long-range 
Coulomb forces between electrons and atomic nuclei exactly, but treats the electrons 
as classical point charges and replaces the Laplacian in the electronic Schrodinger 
equation by a hard-core constraint. Our main results are 

• a simple proof of exact screening in a large- N continuum limit of this model, 
via explicit determination of the minimizer (see (TT5T) ) 

• a proof of approximate screening in this model for large finite N. This is 
done by establishing, via Gamma-convergence, that the discrete minimizers 
converge to the minimizer of the continuum limit (see Theorems 1 1 . 1 \ [TT21 . 

Our model is variational, and describes a system of N particles in M 3 of charge — 1 
('electrons'), with variable positions x±, ...,xn, which Coulomb-repel each other and 
are Coulomb-attracted to M particles of charges +Z a at fixed positions R a ('atomic 
nuclei'): 



Minimize 



N 



1 



(N e N) 



(2) 



i=l l<i<j<N 
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where 



M 

V{X) = 

a=l 



J2 \ x _r\ (Z = (Z 1 ,..,Z M ),Z a >0,R a eR 3 ), (3) 



over the set 

A N := {(x u ..,x N ) G M. 3N \x t - R a \ > d for alH, a} (d > 0). (4) 

The hard core assumption (T4]) may be viewed as a crude "uncertainty principle" 
which prevents electrons from falling into the nucleus, with the hard core radius 
d playing the role of fi. More precisely, the model ([2]), Q, (jl]) arises from the 
full quantum mechanical (Born-Oppenheimer-)Hamiltonian of the electrons in a 
molecule, 



N 

] 



1 N 1 1 

-A + y^ = ^(--A 2i +^))+ J2 T^-y (5) 



X<i<j<N 11 31 



by replacing the one-body operator — hA Xi +v(xi) by the effective potential v e //(^i) := 
f (xj) when \xi — R a \ > d for all i and all a, +oo otherwise. 

The physics of the model (1), (2), (3) is independent of the choice of hard core 
radius d, as long as the hard cores are not overlapping, i.e. 

f(0,oo) ifM=l 
a fc \ (0, lining \R a - R p \) if M > 2. l ° J 

A different choice just corresponds to an overall scale factor of length and energy. 

We proceed to describe our results on the model (2), (3), (4). Comparisons 
with what is known (or expected) in quantum mechanics and related models are 
postponed to the end of this Introduction. 

Our first two observations, the second of which is at first sight somewhat sur- 
prising, are the following: 

Proposition 1.1 a) (Attainment for neutral molecules and singly-negative ions) 



Let Z := J2a=i ^a- F° r N < Z + 1, there exists a minimizer oJVn,z_ on A 
b) (Absorption principle) Every minimizer (xi, ..,Xn) of Vn,z_ on An satisfies Xi G 
Ua=i $a f or a M h where S a = {x G M 3 \ \x — R a \ = d} denotes the sphere of radius 
d centred at R a . 

Proof a) follows from standard arguments in the calculus of variations. The fact 
that unlike in quantum mechanics, attainment can also be shown for N = Z + 1 
comes from the fact that the joint potential exerted by Z particles at x±, ..,xz and 
the nucleus onto an additional particle on a sphere of radius R > max|a;j| is zero 
on average but nonconstant, and hence negative somewhere. We omit the details, 
b) follows from observing that the potential Vn,z[%\-, ■■■,%n) is a harmonic function 
with respect to each particle position Xi and applying the maximum principle. See 
e.g. [Lan72] for a related observation for purely repulsive Coulomb particles confined 
to a bounded (instead of unbounded) set. 
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b) allows to allocate each electron unambiguously to one atomic nucleus. Investi- 
gating the extent of screening means investigating 

- how closely the number of electrons going to any particular nucleus matches the 
nuclear charge (a perfect match corresponds to a zero net monopole moment of the 
atom) 

- how uniformly and symmetrically the electrons distribute themselves around the 
nucleus (this determines the higher net multipole moments). 

In the asymptotic limit when the nuclear charges are large, these questions have 
simple answers. 

Theorem 1.1 Assume d satisfies (EJ), and denote Z := XIqLi Let {x^ , .., x^~^ 
be any minimizer oJVn,z- I n the limit 

N = Z -> oo, ~y -»• z a , (7) 

a) (Neutrality) 

Z a 

b) (Equidistribution) For any measurable Q C S a with area(dQ) = 0, 

N a area(S a ) 

c) (Limit energy) 

Vn,z{x[ N '-\ ..,x { ^- ] ) ly^^l^ \ ^ z a zp 

N 2 2^ d .\R„-Ra\' 



a=l l<a</3<M 

Physically, the results in a), b) and c) are "screening results" which mean, respec- 
tively, that in the above limit 

- the monopole moment of each atom vanishes 

- the higher multipole moments of each atom vanish 

- the interaction energy between the atoms vanishes. 

To understand this interpretation of c), consider, instead of the electronic energy 
Vtf t z_, the total classical energy of the molecule which includes the Coulomb repulsion 
between the nuclei, 



y y 

class • r t/ i \rnuc \rnuc \ v /J a /J ji 

N ,z=™tV N ,z + Vz, Vz = 2. \ Ra - Rp \ 



(Here and below we use the convention that V^ uc = when M = 1.) The formula 
in c) then says that 

rpclass -i M 9 



rpclass M 2 

^N,Z -1- V— > 

~W 2 ' ~2 ^ d 



0=1 

(This is because by (|7j), Z a Zp/N 2 converges to z a zp and hence the internuclear repul- 
sion term cancels the second term appearing in c).) In other words, the limit energy 
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of the molecule just equals the sum of the limit energies of the individual atoms. 
In particular, it is independent of the atomic positions R\,...,Rm, i-e. contains no 
interaction terms. 

The above screening results are a corollary of the following more general result, 
which in addition uncovers interesting behaviour of excess charges moving off to 
infinity in case of negative ions N > Z. 

To include this case we consider, instead of N = Z — > oo, the more general limit 

M N Z a 

N -> oo, Z = Z a -> oo, — -> A, -y- -> z a , (9) 



a=l 



where A G (0, oo) is a filling factor. Positive ions correspond to A < 1, neutral 
molecules to A = 1, and negative ions to A > 1. 

For negative ions, the minimum of Vn,z on ^4.^ is typically not attained (see 
below) and so one needs to relax the restriction to exact minimizers in Theorem 1 1.1 1 
Instead one considers more general low-energy states, in the sense of 

energy difference from infimum << total energy, 

as made precise by the following 

Definition: A sequence {(x^ ,..,x^ )} is called a sequence of approximate 
minimizers of Vn,z_ in the limit ([9]) if 

Vn,z(x[ N '-\ .., x ( ^- ] ) - ini AN V N: z ^ Q ^ 



Theorem 1.2 (Variatonal principle for the limit distribution) For any sequence 
{(x^'~\ .., xff )} of approximate minimizers ofVN,z, in the limit (TJ|) the associated 
measures 

„W3 := Ifx^, (11) 

1=1 

satisfy 

^ fJ'X (12) 

and 

v "- dXi z r N 1 — hM, (is) 

where : .M+(IR 3 \fi) — > ffi U {+00} (see below for notation) is the continuum 
energy functional 



- / K 8\ n E a =i i^u rf M^) + i / V w vh\ d ^ x ) My) if Jd^< A, 

+00 otherwise, 

(14) 

an<i [x\ is its unique minimizer. 
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Existence of a unique minimizer of I\ ;Z is proved in Proposition 2.1 below. Here 
and below our notation is as follows: the halfarrow ^* denotes weak* convergence 
in the space .M(R 3 \fi) of Radon measures on R 3 \ttQ M+(R 3 \tt) denotes the set 
{/i G .M(K 3 \fi) I fi > 0, Q is the union of the hard cores of the nuclei, i.e. Q = 
U^ =1 B d (R a ), B d (R a ) = {x G M 3 | \x - R a \ < d} and z = (z u .., z M ). 

The point about Theorem 1.2 is that the electrostatic continuum energy Ix.z 
which appears in the limit is much simpler than the intricate particle energy V/v,z- 
For neutral molecules (A = 1) or negative ions (A > 1), and nonoverlapping hard 
cores (i.e. flSD); the minimizer of the continuum energy can be determined explicitly, 

a=l 

where H 2 \s a denotes two-dimensional Hausdorff measure restricted to the sphere 
S a = {x G 1R 3 | \x — R a \ — d} (see Proposition 2.1 e) below). This together with 
(fl2|) readily implies the screening results in Theorem 1.1 (see Section 2). 

Theorem 1.2 together with formula (fl5|) also lead to interesting conclusions about 
instability of negative ions. Note that for negative ions (A > 1) the limit measure 
has less mass than the approximating measures, 



J dfi x = 1< A = lim ^ = lim J d^ N ^ ] . 



(The first equality is due to the fact that Yla=i z a = li m Eali if = !•) Physically 
this means that only Z + o(Z) particles stay bound and N — (Z + o(Z)) particles 
move off to infinity. For a precise formulation (as a nonattainment theorem for Vn,z 
when Z is sufficiently large and iV exceeds Z by a nonzero fraction) see Section 5. 

We establish Theorem 11.21 by showing that the particle energy Vn,z an d the 
continuum energy I\ z are related in the mathematically rigorous sense of Gamma- 
convergence, introduced by De Giorgi (see [DM88, Br02] or the beginning of Sec- 
tion 5). Starting point is the observation that the particle energy Vn,z_ can be 
re-interpreted in a natural way as an energy functional on the space A^+(M 3 \fi) 
of nonnegative Radon measures on IR 3 \f2 (with Q as defined below Theorem 1.2). 
Define 



M 

r J \ z i iff i 

V { n,z){v) ■= - z2^F'\ ^-l d v( x ) + o / , I Mx)dfi(y) 

Je. 3 \vi ~y Zj \x — K a \ i J ,/(R3\n)2\diag \ x ~ U\ 



(16) 

if /i = 4X/i=i^i f° r some distinct Xi,..,xn G M 3 \f2, and set V(N,z_){f) '■= +o° 
otherwise. Here and below diag denotes the diagonal {(x, x) | x G M 3 \fi}. Then for 
fi as in the first alternative, we have the identity 



V(Ng)(p>) = -7pV N! z{xi,...,x N ). (17) 



We then show: 



1 Recall that for any closed subset A C R d , M(A) is the dual of the space Co(A) — {/ : A — > 
K | / continuous, f(x) — > for |x| — > oo}, and that a sequence of Radon measures [l v is said to 
converge weak* to fx, notation: \i v — fj,, if J A f d^u — > f A fd[i for all / G C'o(A). 
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Theorem 1.3 (Gamma-convergence) In the limit (fPJ), the sequence of junctionals 
V {N ,z) ■ M + (M. 3 \n) -> MU{+oo} Gamma -converges (with respect to weak* conver- 
gence of Radon measures) to the functional I a,z defined in (T$) . 

Physically, Theorem 11.31 means that the limit functional not just correctly captures 
the ground state energy, but also any energy change by a nonvanishing fraction of 
the ground state energy when the ground state is deformed. 

Note that the restriction to discrete measures has disappeared, and in the domain 
of integration of the second term the diagonal is now included. The latter is essential, 
for otherwise the functional would promote clustering rather than equidistribution, 
and the minimizers would, e.g. in the case of atoms (M = 1), be given by XS X , 
where x is any point on the sphere \x — R\\ = d. In particular an unlimited amount 
of electronic charge could be bound by the nucleus. 

We proceed to compare our results to various results in the literature on other 
models with many-body Coulomb interactions. The attainment result of Propo- 
sition II. lb ) continues to hold in quantum mechanics (Zhislin's theorem, see e.g. 
[Fr03] or see the original Russian article [Zh60]), but requires the slightly stronger 
hypothesis iV < Z + 1 which excludes singly-negative ions. The "no shells" result of 
Proposition II. lb ) is false for true atoms (see [BB55] for experimental data showing 
multiple maxima of the radial electron density in Argon), but interestingly, it is also 
false, e.g., for classical Coulomb particles confined to a disc in two dimensions, in 
which case minimizers would extend into the radial direction [EO00]; but it would 
become true again if the interaction was replaced by the Green's function of the 
two-dimensional Laplacian. For a result related to Theorem 11.21 for repulsive clas- 
sical charges confined to a compact set see [Lan72], where it is proved that every 
sequence of empirical measures of minimizers of the particle system contains a sub- 
sequence converging to a minimizer of the relevant continuum limit. We know of no 
analogues, neither classical nor quantum, of the Gamma-convergence result of The- 
orem [L3j The nonattainment result of Corollary 14.11 for Z large and N > Z + o(Z) 
(see the discussion following Theorem II .2p is known to hold analogously in quantum 
mechanics in the special case M = 1 [LSST], the case M > 1 being open. Numerical 
data of [MDH96] for iV-particle configurations with minimal Coulomb repulsion on 
the sphere suggest that the precise attainment threshold N(Z) of our classical model 
with M = 1 equals Z plus a slowly growing function of Z. For an investigation of 
the higher order energy asymptotics of the latter problem see [KS98]. 

Finally we remark that the passage from the particle energy (T5]) to the continuum 
energy ([TBI) is achieved here via a conceptually new viewpoint which should be of 
more general interest. Instead of parametrizing particle configurations (xi,..,xn) 
with respect to some reference configuration (Lagrangian viewpoint), one considers 
the associated empirical measure const Y^Li S Xi which counts how many particles 
are contained in a given spatial region (Eulerian viewpoint). This allows us here to 
pass to the continuum limit for a frame-indifferent particle system with re-labelling 
symmetry, without any a priori assumptions on admissible particle configurations. 
This strategy should be applicable, at least in principle, to other interesting prob- 
lems, such as many-atom systems interacting via Lennard- Jones-type potentials. 
For why the latter problem, despite involving short-range rather than long-range 
interactions, is in fact harder, see the remarks at the beginning of Section 5.2. 
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Our plan of the paper is as follows. In Section 2 we analyze the limit theory. 
In Section 3 we show how the asymptotic results in theorems 11.11 and 11.21 follow 
from Gamma-convergence and the explicit solution of the limit theory. Section 4 
is devoted to proving nonattainment for sufficiently negative ions, and in the final 
section we establish the Gamma-convergence result of Theorem 11.31 



2 Analysis of the limit theory 

Here we analyze the limit theory (postponing its rigorous justification to Sections 
4 and 6). Its fundamental advantage over the particle system is that it can be 
minimized explicitly. 

Proposition 2.1 (Analysis of limit theory) For any X > 0, z\,..,zm > ; Ri,..,Rm G 
1R 3 ; d > 0, and with Q as defined below Theorem \1.2L the functional 



Hv) = - 1 y\ i — ^s-Mx) + - 

jR^n^lx - Ra\ 2 




2 \x-y\ 



-dfi(x)dfi(y) 



possesses a unique minimizer on M.[o t \] = {jU G .M(IR 3 \f2) | fx > 0, J d/i < A}. 
Moreover when condition (TJP holds, the following statements are true: 

a) (saturation of mass) J dfx\ = min{A, z}, where z := ^2^=1 z &- 

b) (saturation of minimizer) fi\ = fi z for all A > z. 

c) (saturation of energy) e(A) := min^t [0A] / is continuous, strictly decreasing for 
A G [0, z\, and constant for A G [z, oo). 

d) (Minimizer for atoms) If M = 1, 

,,,H 2 \ S , x . zc(X) c(A) 2 

where c(A) = min{A, z} and S is the sphere of radius d centred at R\. 

e) (Minimizer for molecules) If M > 2, X > z (neutral molecules or negative 
ions), 



M tt2\ m 2 



And 2 ' w vy ^2d ^ \R a - R \ 

a=l a=l l<a<f3<M 1 a pl 



where S a is the sphere of radius d centred at R a . 

Note that I is well-defined on A^+(lR 3 \fi) (the space of nonnegative Radon measures 
of finite mass on IR 3 \f2, as introduced below Theorem ll.2p as an element of MUj+oo}, 
because the negative term — YJ Q f R 3\ Q z a \x — R a \~ x djJL{x) is always finite, due to the 



boundedness of the integrand on the domain of integration. 
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The proof of the proposition relies on developing some basic functional analysis 
for the Coulomb self-energy functional 



J{P) -Iff y^—Mx)dii{y) (18) 
2 J Jr6 \x - y\ 

on Radon measures. For smooth, rapidly decaying measures, this functional reduces 
to the very well studied Dirichlet integral for the associated potentials. But this 
standard setting is insufficient here, as the minimizers themselves are singular mea- 
sures which concentrate on lower-dimensional surfaces (see d) and e)). This reflects 
the fact that we are dealing with a support constraint on the measures, as opposed 
to, say, boundary conditions on the potential. 

Denote by C(R 3 ) the set of nonnegative Radon measures on M. of finite mass for 
which J(fi) is finite. Define an extension of J to measures with both negative and 
positive part, as follows: if \i = fa — fa with fa, fa G C(R 3 ), set 

J(fa-fa) := 77 / / i — - — \dfa(x)dfa(y)- [ [ -, — - — -dfa(x)dfa(y) 
2 J J K 6 \x-y\ J J Ra \x-y\ 

■dfa(x)dfa(y). 



2 Jr6 \x — y\ 

Since the first and last term are finite by assumption, and the integrand in the 
middle term is nonnegative, this is well defined as an element of R U {— oo}. The 
key property of J needed in the proof of the proposition is 

Lemma 2.1 J(fa — fa) > for any n\, \i 2 G C(M 3 ), with equality if and only if 
Hi = fi 2 . 

This result is trivial for smooth, rapidly decaying measures, as well as for nonnegative 
measures. That it should continue to hold for rough measures without a sign is well 
known "folklore" in part of the potential theoretic literature (see e.g. [La75]). Our 
proof, given in an appendix, relies on an approximation lemma which concerns 
the behaviour of the Coulomb energy under mollification of measures, and on a 
generalization of an identity of Mattila [Ma95] (see the appendix). 
The lemma readily yields 

Lemma 2.2 / is strictly convex on C(M 3 \fi) ={/ie _M(M 3 \^) | /x > 0, J(fj) < oo}. 

Proof Because the first term of I is linear and the second term is quadratic, we 
have 



2 V 2 J 4 

for any fii, fi 2 in the above set. The assertion now follows from Lemma [2. II 

Finally we will need the following, much simpler, result, which does not rely on 
Lemma 12.11 

Lemma 2.3 I is (sequentially) weak* lower semicontinuous on C(M. 3 \Q), i.e. if h, 
Hj G C(M 3 \fi) with fij ^* fi, then /(//) < lim inf I (fa) ■ 
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Proof This follows, e.g., from Theorem 11.31 and the general fact that T-limits 
are lower semicontinuous (see [Br02, Proposition 1.28]). To keep this section self- 
contained, we include a direct proof, via a simple truncation argument which replaces 
the discontinuous integrand l/\x — y\ in / by a continuous function. Let 

r(*,V):={ ^ .f, 1 *- 8 ;'- (19) 
I a d\x-y\<a, 

and let I a be the functional obtained by replacing the integrand l/\x — y\ in the 
second term of I by f a (x,y). Then liminfj^oo I(fJLj) > lim inf I a (Pj) > Ia(n), 
due to the trivial inequality l/\x — y\ > f a (x,y) and the convergences pij ^* p, 
and fij ® Hj ^* ji <g> p. To finish the proof it suffices to show that lim a _>o I a (p) = 
1(h). If (p ® /•*) (diag) > then this is true because both sides are equal to +oo; 
if (/x®/x)(diag) = then this follows by monotone convergence, because f a (x,y) is 
monotonically increasing in a and tends to l/\x — y\ for all (x, y) ^ diag, and hence 
for (p<g>p)-a,.e. (x, y). 

Proof of Proposition 12.11 Existence of a minimizer is immediate from Lemma 
12.31 any minimizing sequence fi ' is bounded in *M(lR 3 \fi), since J dp^) < A; thus 
there exists a weak* convergent subsequence, by the Banach-Alaoglu theorem; its 
limit must be a minimizer, by Lemma [2.31 Uniqueness follows from Lemma [2.21 

Next we show e). The idea is to first replace the delta functions generating the 
nuclear potential by spherical charge distributions so as to make their self-interaction 
finite instead of infinite; one can then rewrite the energy by "completing the square" . 
So let 

M H 2\ 



p:= > z. 



And 2 ' 

a=l 

with S a as in Theorem 11.11 By Newton's theorem that the electrostatic potential 
exerted by a radial charge distribution onto a point outside it is the same as that 
exerted by the same amount of charge placed at the centre of the sphere, 

M - r 1 

E— ^M=-/ T^-rMv) forallxGM 3 \a (20) 
\x n a \ JR3 \x y\ 

It follows that for any p G A^+(M 3 \fi), 

I(p) = - I I | — - — -dp(y)dp(x) + - I I | — - — -dp(x)dp(y) 
J Jr6 \x-y\ 2 J J Ra \x - y\ 

= \l I r^-Ap-p)(x)d(p-p)^)- 1 - f [ v ^— ] dp(x)dp(y).(2l) 

z j Jmp \ x ~ y\ A j jrs f — y\ 

Since J dp = z and, by assumption, z < A, the measure p — p is contained in M.[o t \y, 
hence by Lemma [2.11 it is the unique minimizer of I on A4\o t \]. 

It remains to evaluate I(p). Using fTSTl) , denoting Pa = z a (4iid 2 )~ 1 H 2 \s a , and 
again using Newton's theorem, 

= 1 f f dp(x)dp(y) 



2 J J \x — y\ 

_l^p r r d Pa (x)d Pa (y) r r dp a (x)dp p (y) _ 1 z 2 a 



x — y\ ^— ' / / |x — y 2 ^— ' d ^— ' LR 
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This completes the proof of e), and also establishes d) in case A > z. 

Next we show d) when A < z. We assume without loss of generality R\ = 0, 
d — 1, and use a symmetrization argument. Given ji G C, let 



voi(50(3)) y SO(3) 

where is the Haar measure on SO (3) and is the rotated measure /j,r(A) = 
/i(-R _1 A) (R G SO(3), A measurable). By the strict convexity of / and Jensen's 
inequality, 

I(p) < with equality if and only if /i — p,. (22) 

Hence it suffices to show that among radially symmetric measures in M.[o,\], the 
unique minimizer of I is given by the formula in d). 

But for such //, we can rewrite I in terms of the radial measure 

v(r) := n(S r ), 5 r = {i G I 3 \x\ = r}, 

as follows. (Note that // G -M[ ,a] is equivalent to v G -M[ a * aZ = {z/ G -M([l, oo)) | z/ > 
0,Jdu< A}, and /i G M A is equivalent to v G .M™cKaZ = {v G 7W([l,oo)) | i/ > 
0, / <iz/ = A}.) By Newton's theorem, 



d/i(y) 



{y\\y\<\x\} F _ 1/1 V ^[l,|rc 



<Mlvl))-n, 

' \x\ 



and by the fact that the electrostatic potential exerted by a radial charge distribution 
onto a point inside it is constant, 



/ 



'{i/||»|>|x|} F-2/1 ^(1^1,00) 12/1 
Consequently, letting Bi := {i G M 3 \x\ = 1}, 



/ = / min{ i ^(N) 

Jr 3 \Bi \x-y\ J[ lj00 ) \x\ \y\ 

and 

= - / -dv{r) + - I I mini-, ■i}di/(r)di/(r / ) =: I(v). (23) 
i[i,oo) r 2 J y [1)0o)2 r r' 

Hence to complete the proof of d), it suffices to show that 

v = min{A, z}5 1 (24) 
is a minimizer of / on A4 r ^f® 1 . But for v G J\/l™ dial ^ since 1 = J du/X, 

Iff ( z : .,11 

[l,oo) a 



I(u) = - / f— r- - ^— + mini-, — }W(r) du(r'). 
1 ; 2./ / rinn yA Ar Ar' V'/i 1 ; 1 ; 



Now for A < z the integrand is minimized pointwise at (r, r') = (1,1), so v{r) = X5i 
is a minimizer on J\4 T ^ dtal . This completes the proof of d). 
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Next, we establish a). For A > z, this follows from the explicit formulae in d), 
e). Thus it suffices to show J d^i\ = A for A < z. The inequality "<" is trivial. 
To prove ">", suppose f d/j,\ < A. Then the measure fl — \i + e(47iR 2 )~ 1 H 2 \ SR (o), 
<3r(0) = {x G 1R 3 | |x| = R}, lies in M.[o,x] for sufficiently small e > 0, and as R — > oo, 
e — > 0, by multipole expansion 

which is smaller than I(fi) for R sufficiently large and e sufficiently small, since 
f dfi — z < \ — z < 0. This contradicts the minimality of /x, completing the proof 
of a). 

b) is a straightforward consequence of a): if A > z and fi\ minimizes I on M[o,\], 
then by a) fi\ G M[o, z ], and hence minimizes / on Ai\o tZ ]. Uniqueness now implies 
A* a = At 2 - 

It remains to prove c). That e(A) is monotonically nonincreasing in A is trivial. 
Continuity at A follows by using the measure \i = (A/A')//y, A' > A, as trial function 
in the variational principle 'Minimize / on .M[o,a]' an d letting A' — > A. That e(A) 
is constant for A G [z, oo) is obvious from b). Finally, strict monotonicity on [0, z] 
can be seen as follows. Let < A < A' < z. By monotonicity, e(A) > e(A'). But 
"=" is impossible, since otherwise uniqueness of minimizers in .M[o,a'] would enforce 
ji\ = fA\r, contradicting the fact that by a), J dp,\ ^ J dfiy. 

The proof of the proposition is complete. 

While from the point of view of Gamma-convergence it is natural to work in Propo- 
sition [2ZD above with the relaxed constrained < J d\i < A (see Theorem IX .3j) . it 
is also of interest to consider the sharp constraint j dfi = A, for this yields a strik- 
ingly simple continuum version of the attainment/nonattainment transition of the 
particle system described in Proposition II. la ) and Corollary 14.11 

Corollary 2.1 Assume condition |6j) holds. Let M\ := {/j G 7W(lR 3 \fi)| / u > 
0, j dfi = A}. If A < z, then the unique minimizer of I on M.\ is given by fi\. If 
A > z, then the infimum of I on A4\ is not attained, the value of the infimum equals 
min_A4 z /, and any minimizing sequence fj,V' converges weak* but not strongly to fx z ; 
in particular J d/i z = z < A = lim^oo J dpP) . 

The proof is straightforward from Proposition 12.11 and Lemma 12.31 See [Ba88] 
for a unifying mathematical setting which subsumes many examples of such at- 
tainment/nonattainment transitions associated with a "loss of the constraint" phe- 
nomenon. 

3 Neutrality and equidistribution 

Here we show how the explicit solution of the limit theory derived above, combined 
with the abstract Gamma-convergence result of Theorem 11.31 allows infer the ab- 
stract convergence result of Theorem 11.21 and the neutrality and equidistribution 
results for the particle system stated in Theorem II. 11 

The main point is that in the topology in which the Gamma-convergence occurs, 
the sequence of (associated measures of) approximate minimizers of the particle 
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energy Vn,z_ is compact. The rest of the argumentation is standard in Gamma- 
convergence. 

Proof of Theorem 11.21 First, note that by (IT7|) and Theorem 11.31 

= inf V (N ,z) - inf / (25) 

in the limit (Q. 

Now let (x[ N ' Z \ ...,x^' Z ^) be a sequence of approximate minimizers of Vn,z- It 
follows from definition (fTUj) and (J25]) that 

Vcat.z)^,-,^)- nif I. (26) 

In addition the associated sequence of measures fj,( N '^ defined in f fTTj) is bounded 
in .M(R 3 \fi) (because ^ N ^ ] > and J RS ^ Q d^ N '^ = N/Z is bounded), and hence 
weak* compact, by the Banach-Alaoglu theorem. 

By the lower bound property (ii) contained in the Gamma-convergence result of 
Theorem 11.31 for every weak* convergent subsequence yiS Nj '—^ the limit fx satisfies 
< liminf I{ir Nh —i'). But by (126]) and Proposition 12. 1| jl must equal the unique 
minimizer of / on -M[o,a]- 

Finally, since every subsequence of fj,( N '—^ converges to this minimizer, so must 
the whole sequence. 

This establishes the theorem. 

Proof of Theorem 11.11 This is a straightforward consequence of Theorem 11.21 
Specializing to iV = Z (whence A = 1) and using Proposition 12. II and the fact that 

= 1 shows that » {N '- ] ^* Eali^ 47 ^ 2 )™ 1 ^ 2 !^ in M(R 3 \n). But by 
Proposition 11.11 the measures fi^ N '—' are supported on U a S a , so the above conver- 
gence also occurs in M.(U a S a ). We now use the well known fact that if a sequence 
of Radon measures fj,^' on any compact (^-dimensional manifold X converges weak* 
to /i, then for all Borel sets A C X with n(dA) = 0, ijP\A) — > [i{A). Consequently 
for A C S a 

§{ x i N) I x< i N) e A} _ f xAxW N,z ){x) ^ f xAx)Mx) 



N J Sa AX/ W Js a AW W area(S, 

This proves Theorem 11.11 



4 Instability of asymptotically negative ions 

Here we show how Theorem 11.21 together with the saturation of mass phenomenon 
of Proposition 12.11 a) leads to a nonattainment result on Vn,z which complements 
Proposition II. lb ) . 

Corollary 4.1 (Instability of asymptotically negative ions) 

Let N*(Z_) = supjiV e N| inf a n V^,z attained}. Then as Z = z^2 a=1 Z a ~^ °°> 
Z a / Z — > z a , 

Z 
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Proof The simple attainment result in Proposition 11.11 a) implies that N*(Z_) > 
Z + 1; hence it is clear that liminf^^oo N *^r* > 1. The nontrivial assertion in the 
corollary is that 

NJZ) , . 

limsup K ~ J < 1. (27) 

Z^oo Z 

But this is a direct consequence of Theorem 11.21 and Proposition 12.1a ). as follows. 
Denote the value of the limsup in (I2TI) by A*, and consider any subsequence real- 
izing it, i.e. N m (ZV)/ZV> -> A* as j -> oo. Abbreviate iVj := iV*(Z (i ), and let 
(x^, ...,x^) be a minimizer of V^. z u)- 

On the one hand, by Theorem 11.21 Proposition 12.1a ). and the fact that Yl a z a = 
1, the associated measure satisfies .M(IR 3 \fi), for some measure \x 

with J k3 , n <i/i = min{A*, 1}. 

On the other hand, by the absorption principle (Proposition 11.11 b)), letting 
S := D a S a , J s dfx^'— 3 ^ — > A* and, due to the fact that the above weak* convergence 
also occurs in A4(S), J s dfi = A*. 

Consequently minjA*, 1} = A*, or equivalently A* < 1, as was to be shown. 

5 The continuum theory as Gamma limit of the 
many-particle Coulomb system 

Here we show that the continuum theory (|T4"|) arises in a mathematically rigorous 
way (namely as a Gamma-limit) from the many-particle Coulomb system, i.e. we 
prove Theorem 11.31 As emphasized in the Introduction, despite a result of this kind 
being — in our view — very natural, we know of no case of any many-atom or 
many-electron system where such a result has been previously established. 

Recall (e.g. from [Br02]) that a sequence 1^ :X-^lU {oo} of functionals on 
a topological space X is said to T-converge to J : X — > IR U {oo} if for all /i G X 
we have: 

(i) (Ansatz-free lower bound) For every sequence \ij G X converging to \i we have 

< liminfj_ 1 . 00 /^ ) (//j). 

(ii) (Attainment of lower bound) There exists a sequence //j G X converging to \i 
such that = limj-^oo I^'(fMj). 

In our case, X = M.[o,x] (see Theorem II .31) . the space of nonnegative Radon measures 
on M 3 \f2 of mass < A, endowed with the weak* topology, and = /( Ar(3) >^ (j) ) , where 
N U) G Nj z® = ( Z C0 Z $), Z U) = J2™ =1 Zg\ and 

zi j) 

-^-y -> A G (0, oo), -^-y -> 2 a (j -> oo). (28) 

In order to establish Theorem 11.31 we need to verify (i) and (ii). 
5.1 Proof of the lower bound (i) 

Suppose that Hj ^* fi. We may assume without loss of generality that I^\fXj) < oo 
for all j (because if J := {j G N | I^'ijij) < oo} is finite the assertion is trivial and 
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if it is infinite then timinfj_ foo i'^(/4 7 -) = liminfjej ) j_ +00 /^(/i J -)). By passing to a 
subsequence we may in addition assume that -> liminf^oo J (j) (^)- 

We use a truncation argument as in the proof of Lemma [2731 We let f a , I a be 
as defined there, and compute 



r -\ iff i 

^ (fc) (/ifc) = - / V - p r ^fc(^) + 77 / / , I \dn k (x)dn k (y) 

J^oTi W \ x - R a\ 2 J </(R 3 \c) 2 \diag \ x ~ v\ 

. M ^(k) Iff 

> ~ Y] 7ft) i a p r d(i k (x) + - / / f a (x,y)dn k (x)dn k (y) 
Jm 3 \n ~y Z ( >\ x ~ tta\ 2 J J(R3\n)2\diag 

M z { k ) . 



A/ z {k \ f 1 

Tdfi k (x) 



\x-R a \ ^ K ' 2(Z( k )) 2 a' 

Using the fact that the last term on the right hand side and the factors in the middle 
term tend to zero by (jUJ) and that if ^* /x in .M(IR 3 \f2) then /x fc <g> ji k ^* /x (g> /x 
in M. ((M 3 \f2) 2 ), letting tend to infinity gives 

limlW(/x fc ) > lim J Q (/x fc ) = I a (/x). (29) 

But as shown at the end of the proof of Lemma 12.31 lim Q _>o / a (/x) = This 
establishes (i). 



5.2 Proof of the upper bound (ii) 

Fix a sequence iY» -> 00, satisfying Given /x G A^ + (R 3 \Q), we need to 

construct a sequence /x„ G A^+(IR 3 \f2) (or, in Gamma-convergence terminology, a 
"recovery sequence") such that /x„ ^* /x and limsup^^ l^ N{v) ^ M )(/x„) < /(/x). 

This is achieved by a careful multiscale construction, by introducing a mesoscale 
h with particle spacing << mesoscale << diameter of support of /x (see Step 2) 
and approximating /x in each mesoscale region by a suitable number of Dirac masses 
placed on some suitable lattice (see Step 4). The number of Dirac masses is governed 
by the amount of mass to be accommodated in the region (see Step 3). 

A first difficulty is that unlike in usual arguments establishing density of discrete 
measures, the amplitude of each Dirac mass is fixed exactly to be 1/Z, which leads 
to a mass error of order 1/Z in a typical mesoscale region. 

A second, and more fundamental, difficulty is that one expects the energy to be 
highly sensitive to the precise placement of the particles; but the precise structure 
of approximate or exact minimizers of the many-body Coulomb interaction is un- 
known mathematically. We know of no attempt to prove minimizers are cystalline, 
or approximately crystalline, let alone to establish the optimal lattice structure - 
neither for Coulomb interactions nor for any other realistic interaction law in three 
dimensions. 

Very remarkably, the long range nature of the Coulomb force, usually considered 
a complicating rather than a simplifying feature, works in our favour. It implies 
that the energy is dominated by long range contributions, and so at short range a 
rough knowledge of bondlengths (to within a factor) turns out to be sufficient, as 
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long as knowledge of the long range distances, governed by the "packing density" , 
is precise. The key point in the proof is the implementation of these ideas in Step 
5. 

Step 1 Reduction to compactly supported measures of finite energy, mass X, and 
bounded Lebesgue density 

If J(/i) = oo existence of a recovery sequence is trivial: for instance the sequence 
jj, u = ji will do. So we may assume < oo. 

By a standard approximation argument, we may also assume that fi has compact 
support in R 3 \f2. 

A little less trivially, we claim that it is enough to establish existence of a recovery 
sequence for measures with J dp, — A. This is because if f dji > A, then I(fi) = oo 
and we are back in the case dealt with above, whereas if J dfi < A we can always 
"place unwanted mass at infinity". More precisely, if J dfi = A < A, we choose 
N v € N, N v < N v , such that N v /Z v — > A, apply existence of a recovery sequence fi v 
for fi with respect to the functionals I^ Nv ' Zv \ and set fi u = fi u + '^2k=x N " {^/ Z u )5 x k 
with minfc \x^\ — > oo, min^ \x k v — x v \ — > oo. It follows that fi v — fi v ^* and 

Finally we claim that we may assume that \i has bounded Lebesgue density, 
i.e. dfi(x) = m(x)dx for some m G L°°(R 3 \fi). (We thank Stefan Miiller for this 
idea, which facilitates the simple energy error estimate (jUJ) via the uniform bound 
(J40]) on the local lattice spacing below. It replaces our original more complicated 
energy estimate via an integral bound on the lattice spacing.) This is because given 
any measure /i G .M + (R 3 \fi of mass A and compact support, one can construct a 
sequence of measures /i^ G .M + (R 3 \fi) fl L°°(R 3 ) of mass A and compact support 
with and liminf e ^o Ht 1 ^) — Ht 1 )- Indeed, we claim that the mollified 

measure 

„«■>:=*..„=/ #.(,)„(-*)* 



with <f> e (z) = e" 3 0(e- 1 2), <p G Cg°(5i(0)), <p > 0, / R8 <p = 1, has the required 
properties. For the elementary proof of weak* convergence to \x see e.g. [M] Theorem 
1.26. To verify that liminf e ^ I(^) — ^(/-Oj we write /(//) = f R3 ^ n v(x)dfi(x) + 
J(/i), with J as in (jTSl) . The weak* convergence of //^ implies J" R3 , Q v(x)d(j,^(x) — > 
/r 3 \q v ( x )dfi(x), and the translation invariance and convexity of J and Jensen's 
inequality imply 

J(n)= f (f> e {z)J{n{- - z))dz> j( I <l> e (z)fi(- - z)dz) = J(/i (e) ) 

for all e, establishing the assertion. We remark that by weak* lower semi-continuity 
of J, in fact one has I(fi^) — > but this is not needed here. 

Step 2 Discretization o/R 3 \f2 into mesh of size h 

From now on we fix a measure fi G .M+(R 3 \f2) which has finite energy and is 
compactly supported in R 3 \f2. Hence for ho sufficiently small, supp /i is contained 
in some finite union of disjoint cubes of sidelength ho, Q' = U^ 1 Q/ l0 (i?^,), with 
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Q' C R 3 \fi, where Qh{R) denotes a cube centred at R of sidelength h, {x G M 3 \ — 
h/2 < (x — R) ■ tj < h/2}, and where ei, e-i y are the standard basis vectors of M 3 . 

Now given any ngN (to be chosen later, depending on N and Z) we obtain a 
mesh of size h := h /n by dividing each cube of sidelength /i m to n 3 smaller cubes 
of sidelength ho/n. This way we obtain a disjoint family of cubes {Q® Q f 
sidelength h whose union contains supp \i. 

Step 3 Choice of number of Dime masses in each region and mass error analysis 
For given /i, N, Z, we need to approximate //|n(o by a measure of form 



(N,Z) 



Lj G N U {0}, ^TL, 



AT, 



(30) 



fc=i 



because otherwise ^(jv,z)(/-* ~ ) is infinite. In particular, the allowed mass has to 
be an integer multiple of 1/Z, enforcing a mass error. Here we deal with the latter, 
postponing the choice of positions x\ k ^ to the next step. 

It will be convenient to introduce a trivial amplitude factor 



(31) 



(note (f)(N, Z) — > 1 in the limit ([2])) and approximate the measure • /i, because 

Choose £; G N U {0} such that is a good approximation to the mass of (ft/j, 
in gW, i.e. 



< 



+ i 



(32) 



This implies EA/Z < N/Z < E^qM)^ + l )l Z - Hence (J2A) + r = N for 
some integer r less or equal to the number of indices i with fj,(Q®) > 0. It follows 
that if for r such indices we set Li :— t{ + 1, and let Lj := otherwise, we have 
YU Li = N,Li = when MQ (i) ) = 0, and, by 



; " 1 < MQ W ) < 7 1 



(33) 



Z - rr ^ '- ^ 
Hence if ^ N, ^\n{i) is given by f[3"Ul) . regardless of the choice of positions xf 1 ' we have 

\A\ 



Ad((j>fi) 



< — for all A G 

z 



(34) 



Let Cj denote the centre of the cube QW. Then by (134"1) . if g G C(M 3 \f2) and 5 is its 
modulus of continuity on the lengthscale of the mesh, 



5 := sup 

\x-y\<maxi diaULQ« 



g(x) 



g{y) 



then 



g{(f)dii) 



R :i \n 



gdfi^ 
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\J2 j ^g{x)-g{ Ci ))d{^){x) + g{ Ci )(d{^){x)-d^ N ^\x)) + (g{ Ci )-g{x))d^(x)) 
<sfd{^) + + sfd^ = 24 + supM§| (35) 



Z " iai h 3 Z 

where the factor N^h^/h 3 = N^n 3 in the last term is the number of cells Q^. 

If N — > oo, Z — ► oo, /i = /i(iV, Z), it follows that the right hand side tends to 
zero for all g G Co(IR 3 \f2) provided the meshsize h satisfies 

h^O ("small mesh") (36) 

(whence 5 — > 0), 

£-1/3 

— r y ( "particle spacing smaller than mesh" ) (37) 

(whence the second term in (155]) tends to zero). Hence if ( 136]) and (157}) hold, then 
^{n,z) _ ^ ^* o an d hence, thanks to 0(JV, Z) -> 1, ^* \i. 

To understand the meaning of (157]) . it is instructive to consider the case when /i is 
the uniform measure on some region of finite diameter, and when l£ N & is positioned 
on a periodic lattice in this region. Because jjS- N '^-> has iV Dirac masses, the lattice 
spacing must be ~ N~ l l 3 ~ Z~ l l 3 and hence Z~ x l 3 /h ~ (particle spacing) /meshsize. 

Step 4 Choice of positions of Dirac masses 

The Dirac masses in each will be positioned by placing a lattice of spacing 
~ (maxj Lj) -1 / 3 in QW. First, we estimate max^Lj, using the fact that jx has 
bounded Lebesgue density. By (133]) and the fact that li(Q^) < h 3 \ 

L,<Z0(iV, Z)||/i|| L o / i 3 + l. 

Hence since h 3 Z — > oo and <p(N, Z) — > 1, we may assume 

L t <Ch 3 Z =: L+, C r = 2||/i|| i ac. (38) 

Now we choose a lattice aZ 3 of sufficiently small spacing a so that each contains 
at least Lj lattice points. Because has sidelength h, it suffices to take 

h , \ 

(39) 



r(L + )V3] > 



where [(L+) 1 / 3 ] denote the smallest integer > (L+) 1 / 3 . For future use we note that 
since L + > 1, \(L + )^ 3 ] < 2{L + ) 1 ' 3 and so 

1 h hi 

<a<—jr= ^ 1/a „ 1/a - (40) 



2C 1 /3^i/3 2Li /3 ~ ~ L+ /3 C X I 3 Z X I 3 ' 

Now for each QW, we choose {xf~\ a;^ } as a subset of aZ 3 H Q^, and let fi^ N, ~^ 
be as defined in (1301) . 

Step 5 Analysis of energy error 
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With the above choice of the xf\ with h = h v chosen to satisfy (1361) and (1371) . and 
with defined by (1391) . (!38|) . we claim that 

lim/ (^,^) (/1 (*&) ) = /(/i) . (41) 

By the weak* convergence of y^ Nu ^v) to /x, we can immediately pass to the limit in 
the electron-nuclei interaction: 

v{x)d[i {JS!v, - v \x) — > / v(x)dfx(x). 



As for the electron-electron interaction, we decompose 



2 J J(K 3 \n) 2 \diag F ~~ 2/1 2 J J(iR3\n)2\diag I s — 2/ 

' f f fJd^ N ^®d^ N »^-d l i®d^ ' ;V " 



2 7 7(r3\q)2 V ' / 2a (Z u ) 2 

+U [ (r^—i - /«W"^ ® dfJL^ -Uf (-^— - f a )d» ® dp. 

2 7 i(R3\ f7) 2\diag v F -2/1 y 2 y Jgg3\Qy \ \x — y\ J 

Because the measures ^ Nv '—^> are supported in Q! , and / Q is continuous on the 
closure of Q! x Q! , it follows from the weak* convergence of ^ Nv '—^ that the first 
term tends to zero as v — > oo. The second and fourth term are < 0, and the 
integrand in the third term is bounded from above by \^_ y \ X{\x-y\<a\- It follows that 



limsup(/ { ^'^ ) (/^ ( ' V, "- J )-^(^)) <limsup~ J ^" R3vn)2vdiag ^—dfi^'^dfi^ 

n{\x-y\<a} 

(42) 

for all a > 0. 

It remains to show that the right hand side of (l4"2l) tends to zero as a — > 0. By 
the trivial estimate for any nonnegative Radon measure m 

1 i m ® im <(l dm \ • sup / r^dmfo) 



( R 3\n)2\diag |x - y 
and the fact that J d^ Nv '—^ — > A, it suffices to show that 

limsup sup / 1 d^ N »^\y) -» (a -»• 0). (43) 

z/^oo xesupp AtC^"'^) JyeB a (x)\{x] F ~ V\ 

We estimate the right hand side by neglecting the fact that only a subset of a v l? 
carries Dirac masses and applying the lower bound (1401) on a v : 

sup / ^— d ^\y) < z; 1 sup Yl r^- 



xesupp^-^) JyeB a (x)\{x} F - 2/1 !BeavZ3 »€e v z8\{-}' a: ~ 2/1 

Pl{ I y — a; I < q } 

< 8C ■ £ -1. (44) 



S/Ga,Z3\{0}, |i/|<a 
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2/1 



Now the latter is a Riemann sum and since its "meshsize" a v < (CZ U )-^ 3 -> by 
( BUI) , it converges in the limit v — > oo to 8C f , <Q \y\~ 1 dy. But explicit evaluation 
of this last integral shows that it tends to zero as a — > 0. This establishes ( B3l . 
and completes the proof that the constructed multiscale lattice measures ^ Nl "~ v 
constitute a recovery sequence. 

The proof of the Gamma-convergence result Theorem 11.31 is complete. 

Appendix: The Coulomb norm on Radon measures without 
a sign 

We show here that the Coulomb energy on Radon measures without a sign is strictly 
positive definite, and hence gives rise to a norm. This was used to determine the 
minimizer of the continuum energy in Proposition 2.1. The result is standard "folk- 
lore" in potential theory (see e.g.[La72]), trivial for measures which are either smooth 
and rapidly decaying or nonnegative, and surely well known to experts, but we were 
unable to find a reference. 

Our proof is based on two lemmas. The first is an approximation lemma for 
Radon measures concerning the behaviour of the Coulomb norm under mollifica- 
tion; analogous statements are very well known in Sobolev spaces. The second is a 
representation formula via the Fourier transform which is a modest generalization 
of a corresponding identity in [Ma95, Ch.12] for nonnegative Radon measures with 
compact support; our argument is different as the restrictions on sign and support 
are important for the argument in [M a95] . 

Let .M+(R 3 ) denote the set of nonnegative Radon measures of finite mass on R 3 , 
and let C(R 3 ) denote the subset of such measures /i whose Coulomb energy J(/i) 
defined in f fl8|) is finite. Define an extension of J to measures with both negative 
and positive part, as follows: if fi = fii — \ii with jjli, // 2 G C(R 3 ), set 

J(//i-// 2 ) := - / / i — - — rd/j ll (x)diJ ll (y) - [ [ - — - — rd^xjd^iu) 
2 J Jus \x-y\ J J R 6 \x-y\ 

-dfx 2 (x)dfx 2 (y)- 



2 J Jr6 \x — y\ 



Since the first and last term are finite by assumption, and the integrand in the 
middle term is nonnegative, this is well defined as an element of R U { — oo}. Denote 
by C*o(R 3 ) the function space {u : R 3 — >• R | u continuous, u(x) — > as \x\ — > oo}. 

Lemma Al Let G C*o(R 3 ) be nonnegative, radially symmetric, and satisfy j R3 <fi = 
1. Let Hi, /i 2 G C(R 3 ), let \i = — Hi, ^nd for e > let // denote the mollified 
measure 

H e (x) = (<f> e * n)(x) = / <fr € (x - x')dn(x'), (45) 



where <f> e (z) = e 3 (p(e 1 z). Then J (fx) = \im £ ^ J([i e ) ■ 

Proof It suffices to consider the middle term in the definition of J, i.e. to show 

i -dfii{x)dfi2{y) = lim / / -dii\(x)dii £ 2 (y) (46) 

R6 F ~ V\ e ^°J Jm 6 \ x - y\ 
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(for the other terms, set fi\ = fx?). By Fubini's theorem, 



M l (x)d^ 2 (y) = 11(11 MX y,) dxdy)d^x'W 2 {y'). 
J JRa V Jm.6 \x — y\ / 

(47) 

Clearly the term in brackets convergers pointwise to l/\x' — y'\ as e — > 0. Moreover 
setting z = x — x' and using the radial symmetry of <p E together with Newton's 
theorem that the electrostatic potential exerted by a spherical charge distribution 
onto a point inside it is constant, while that exerted onto a point outside it is the 
same as that exerted by the same amount of charge placed at the centre, 

<f) £ {x-x') f (f) e (z) f . 1 1 un , 
-ax = / —dz = / mm{ — , , -\<pAz)dz 



x~y\ " J R 3 \z-(y-x')\ ' J R3 \z\' \y-x 

< I ^{z)dz = -r-^—fr. (48) 

Multiplying by <p E (y — y'), integrating over y and applying (PfB]) again yields 

(j) E (x - x')(f) £ (y - y') /" 1 , 1 
; ; dxdy< / -j -cj) E (y-y)dy < -— -. 

f - y\ Jrz \y - x \ f - y \ 

Hence by dominated convergence, the right hand side of (HTj) converges to the left 
hand side of f|46l) . This establishes the lemma. 

Next, we derive an expression for the Coulomb energy in terms of the Fourier 
transform, defined for any nonnegative Radon measure of finite mass, /i G .M+(R 3 ), 
by 

j2(k) := / e- ik - x dfi(x). 




Note that for any such measure, its Fourier transform is a bounded continuous 
function. 

Lemma A2 Let fi 2 G C(R 3 ). Then 

Air 



J(/il ~ ^ = 2 (2^ L " ^ ?dk - m 

Proof We approximate fi±, fi 2 by smoother measures, as follows. For e > 0, let 
(j) £ (x) := (2n£)~ 3 / 2 e~ x ^ 2e \ Given any fi G C(R 3 ), denote the associated mollified 
measure (HSD by ff. Then // G L^M 3 ) nL°°(R 3 ), because / <p £ *n = (J <j> e )-(J dp) = 
J d\i < oo and sup(0 e * fi) < (sup0 e ) • J dfj, < oo. In particular // G L 2 (IR 3 ), and 
so its Fourier transform fi £ is well-defined as an element of L 2 (IR 3 ). By standard 
Fourier calculus and the fact that f(x) = l/\x\ has Fourier transform 4ir/\k\ 2 , 

f? = 4> £ ■ jl, 

11 /" 4tt ^ 
m = 2(2^/ R3 W mk)?mk) (50) 

As e — > 0, because 4> E {k) = e~ efc2 / 2 we have that \4> E (k)\ 2 tends monotonically to 1, 
and hence by monotone convergence the right hand side of (1501) tends to the right 
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hand side of (jUJ). On the other hand, by Lemma Al the left hand side of (150]) tends 
to the left hand side of ( 1491) . This establishes the lemma. 
Finally we assert: 

Lemma A3 Let fix, ^2 G C(R 3 ). Then J(/ii — ^2) > 0, with equality if and only if 
/ii = /i 2 . 

Proof Nonnegativity is clear from ( 14*91) . Moreover the right hand side of ( 149]) is 
strictly positive unless fli = 1I2 Lebesgue-almost everywhere. But by continuity of 
the fii, this means fli = [2 2 , and hence /ii = /i 2 . The proof is complete. 
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